Scalable quantum search using trapped ions 
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We propose a scalable implementation of Grover's quantum search algorithm in a trapped-ion 
quantum information processor. The system is initialized in an entangled Dicke state by using 
simple adiabatic techniques. The inversion-about-average and the oracle operators take the form 
of single off-resonant laser pulses, addressing, respectively, all and half of the ions in the trap. 
This is made possible by utilizing the physical symmetries of the trapped-ion linear crystal. The 
physical realization of the algorithm represents a dramatic simplification: each logical iteration 
(oracle and inversion about average) requires only two physical interaction steps, in contrast to the 
large number of concatenated gates required by previous approaches. This does not only facilitate 
the implementation, but also increases the overall fidelity of the algorithm. 

PACS numbers: 03.67.Ac, 03.67.Lx, 03.67.Bg, 42.50.Dv 



-I— > 

c ■ 

=5 

cr 



> 

m 

O 
O 



X 



I. INTRODUCTION 



One of the most celebrated applications of quantum 
information processing is Grover's quantum search algo- 
rithm, which allows an initially unknown element to be 
determined from Af equally likely possibilities in 0{\fl\f) 
queries [l| . This outperforms the optimum classical strat- 
egy (a random 'trial and error' of elements), which re- 
quires O (Af) steps on average. In addition to providing 
a speed-up of the unstructured search problem, Grover's 
algorithm can also be adapted to look for solutions to 
a range of mathematically difficult problems that have 
structure, by nesting one quantum search inside another 
0. As with other applications of quantum computing, 
the benefits of a quantum over a classical approach in- 
crease with the size of the database. Indeed, it has been 
suggested that the primary resource for quantum com- 
putation is a Hilbert-space dimension, which grows ex- 
ponentially with the available physical resources Q . An- 
other important consideration is that the physical and 
temporal resources required to implement the algorithm 
do not grow too rapidly with the register size Af. 

Proof-of-principle quantum search has been success- 
fully demonstrated in nuclear magnetic resonance [H-[(| , 
linear optical and trapped-ion systems @, as well as 
with individual Rydberg atoms [To[ and in classical optics 
[IH . Of these, only the trapped-ion platform possesses 
a fully scalable Hilbcrt space and in this sense it is real- 
istically the only candidate for performing a practically 
useful quantum search. We note, however, that while 
the trapped-ion system is scalable [L2\ , the largest di- 
mensional quantum search so far performed with trapped 
ions was for a database size of Af = 4 0. Extending the 
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approach of Ref. Q to a large number of ions is highly 
demanding, since it requires the ability to construct, with 
very high fidelity, a great number of multiply-conditional 
gates. For two ions, as in Ref. @, each oracle query 
amounts to a controlled phase gate between the ions. 
When N ions are involved, however, the oracle operator 
becomes a gate which is multiply-conditional upon the 
internal state of all N ions. Although such a multiply- 
conditional gate may be decomposed efficiently using a 
series of one- and two-ion gates [l3j . in practice such a 
synthesis becomes a daunting task, even for a moderate 
number of ions. For example, for a four-qubit register, 
the oracle may be constructed using 13 two-qubit con- 
ditional gates [l3[ , each of these requiring 5 consecutive 
physical interactions on average, which is beyond the ca- 
pabilities of current experiments (l4T4T6j . The above ex- 
ample makes it apparent that there is a clear distinction 
between an approach, which is only formally scalable, and 
an approach, which is realistically scalable using current 
technology. 

It is highly desirable, therefore, to find new ways of per- 
forming quantum search in a scalable system, which docs 
not require an exorbitant number of elementary physical 
interactions. In this paper we propose a novel approach 
to perform Grover's search algorithm. The linear crystal 
of N ions is prepared in a symmetric Dicke state with 
N/2 excitations, \W^ 2 ); hence the register dimension 
scales exponentially with the number of ions. The imple- 
mentation of the quantum search involves a series of red- 
sideband laser pulses addressing all, or half, of the ions in 
the trap. Thus each of the inversion about average and 
the oracle query can be produced in a single interaction 
step. Consequently, the total number of physical steps is 
the same as the number of algorithmic steps. 

The implementation of the Grover's search algorithm 
proposed below follows earlier proposals, wherein the 
database scales linearly [l8| (see also [2(|) and quadrati- 



2 



\w) = -j=y 

1 1 VM^ 



cally [19( with the number of ions N. The database size tion of all basis elements, 
here, which is M ~ 2 N ^2/ (nN), scales exponentially 
with iV and it is therefore of fully quantum nature. The 
reduction of the full Hilbert space of dimension 2 N by 
a factor of \J~N is compensated by the absence of any 
ancilla qubits or the need of error correction because of 
the ultra-compact implementation, in which each logi- 
cal mathematical step is realized with a single physical 
interaction (laser pulse). 

The remainder of this paper is organized as follows. 
A brief review of Grover's algorithm for a quantum me- 
chanical speed-up of the unstructured search problem is 
given in Sec. |H] Particular attention is given to the 
ideas of amplitude amplification and generalized reflec- 
tions, which lie at the heart of the search algorithm. Sec. 
IIIII introduces a model Hamiltonian describing a single 
laser pulse illuminating a chain of trapped ions, and a 
convenient Hilbcrt-spacc factorization, with which to de- 
scribe the dynamics. By making use of the simplified dy- 
namics in this factorized basis, we construct recipes for 
the approximate synthesis of the two reflections required 
for a quantum search, and examples of numerical simula- 
tions are presented in Sec. IIV1 In Sec. |V]wc summarize 
our findings. 
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Central to the operation of the quantum search al- 
gorithm is the idea of generalized complex reflections, 
known in the computer science literature as Householder 
reflections (HR) |21| : 



M^) = l + (e*-lM<# 
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When the phase is set equal to tt, the effect of M^(<f>) 
on any vector is to invert the sign of the component of 
this vector along while leaving all other components 
unchanged, which amounts to a reflection with respect 
to an (J\f— l)-dimensional plane orthogonal to Al- 
lowing <j) to take arbitrary values extends the concept 
of reflection, by imprinting an arbitrary phase onto the 
component along \ip), rather than a simple sign inversion. 
Householder reflections are widely used in classical data 
analysis and also constitute a powerful tool for coherent 
manipulation of quantum systems (22| - [26| . The core com- 
ponent of Grover's algorithm is a pair of coupled House- 
holder reflections, which together form a single Grover 
operator G: 



G = M w (tp w ) M„(tp„). 



(4) 



II. GROVER'S SEARCH ALGORITHM 

Grover's algorithm provides a method for solving the 
unstructured search problem, which can be stated as 
follows: given a collection of database elements x = 
1,2,..., Af, and an oracle function f(x) that acts dif- 
ferently on one marked element s to all others, 
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find the marked element in as few calls to f{x) as possible 
If the database is encoded in a physical system that 
behaves classically, then each oracle query can only act 
on a single database element. In this case, the optimal 
search strategy is simply a random selection of elements; 
on average, it will be necessary to make approximately 
Af/2 calls to the oracle before the marked element s is 
located. The idea underlying Grover's algorithm is to 
encode the database in a physical system that behaves 
quantum mechanically. Therefore, each possible search 
outcome is represented as a basis vector \x) in an TV- 
dimensional Hilbert space; correspondingly, the marked 
element is encoded by a marked state \s). Hence one can 
apply unitary operations (involving the oracle function) 
to superpositions of the different search outcomes. It is 
thus possible to amplify the amplitude of the marked 
state |s) using constructive interference, while attenuat- 
ing all other amplitudes, and locate the marked element 
in 0(VW) steps. Before the execution of the algorithm, 
the quantum register is prepared in an equal superposi- 



According to standard nomenclature, the operator 
M s (ip s ) is referred to as the oracle query, while M-w&w) 
is known as the inversion-about-average operator. 

We note that with the initial state given in Eq. (|2|), 
and during successive applications of the operator G, the 
state vector for the system begins and remains in the 
two-dimensional subspacc defined by the non-orthogonal 
states \s) and \W). Each application of G amplifies 
the marked state population until it reaches a maximum 
value close to unity after uq iterations, at which point 
the search result can be read out by a measurement in 
the computational basis. 

The problem of how to optimize the quantum search 
routine by allow ing arbitrary ipw and (p s has been stud- 
ied extensively [27H29( | . It is found that the maximum 
possible amplitude amplification per step of the marked 
state arises when the phases tpw and tp s are both set to tt 
(as in Grover's original proposal [l[). The corresponding 
minimum number of search steps ng m is given by: 



(5) 
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where [n] denotes the integer part of n. However, this 
choice of phases is not unique. For large Af, as long as 
the phase matching condition ipw = tp s — <p is satisfied 
[28j , a high fidelity search can be performed for any value 
of <p in the range < ip ^ tt and for certain values of ip, 
a deterministic quantum search is possible [27j ■ 

In the following two sections, wc shall describe how to 
perform a physically efficient quantum search using the 
dynamic symmetries in a system of trapped ions. 
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III. ION-TRAP IMPLEMENTATION 

A. Model Hamiltonian 

We consider JV ions confined and laser cooled in a lin- 
ear Paul trap, each with two relevant internal states |0) 
and |1), with respective transition frequency ujq. The lin- 
ear ion crystal interacts with a laser pulse with frequency 
cjl = wo — cJtr — S, where Wt r is the axial trap frequency 
and 5 is the laser detuning from the first red-sideband 
resonance. We assume that the phonon spectrum can 
be resolved sufficiently well that only the center-of-mass 
mode is excited by this interaction and that other vi- 
brational modes can safely be neglected. In the Lamb- 
Dicke limit, and after making the optical and vibrational 
rotatin g-w ave approximations, the interaction Hamilto- 
nian is 1301 : 



JV 



H I (t) = hJ2 



k=l 



g(t) {(ria + a k a 1 ) + -a k 



(0) 



Here a k = |lfc)(0fe| and a k = |0fc)(lfc| are the raising and 
lowering operators for the internal states of the fcth ion, 

a k is the Pauli spin matrix and d T and a are respectively 
the creation and annihilation operators of center-of-mass 
phonons. The coupling between the internal and mo- 
tional degrees of freedom is g(t) = rjCl(t)/2y/N, where 
i] = -\/fi.|k| 2 /2Ma;tr is the single-ion Lamb-Dicke param- 
eter, with k being the laser wave vector and M is the 
mass of the ion. The function f2(t) is the real-valued 
time-dependent Rabi frequency. 

The Hamiltonian (j6|) conserves the total number of ex- 
citations (rii ionic plus n p motional) , which in the scheme 
we propose, is half the number of ions, i.e. N/2 (with iV 
even), rii + n p = N/2. The energy pattern splits into 
manifolds corresponding to rii ionic and n p = N/2 — rii 
motional excitations. Each manifold is C^-fold degener- 
ate, where C N = N\/rii \ (N — m)l. It is readily verified 
that the dimension of the rii = N/2 manifold, T>, grows 
exponentially with JV; indeed, for large JV we have 
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+ C(JV" 2 ) 



(7) 



The subspace of the overall Hilbcrt space, which spans 
the manifold T>, we shall use to represent the state of the 
register in Sec. IIV) T> is the set of states, which encode 
the database of N elements. 



B. Hilbert space factorization 

It is possible to move to a new basis in which the 
Hilbert space is factorized into a collection of noninter- 
acting chains of states. The new states we shall call 
MS states, since they can be obtained by the multilevel 
Morris- Shore (MS) factorization [11]. 
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FIG. 1: The MS basis states for N = 6 ions with 3 excitations 
consists of a series of independent chainwise linkages. Since 
the total number of excitations is half the number of ions, 
then from Eq. ((9]) follows that states with mj > are inac- 
cessible and are therefore not shown. The states that make 
up the longest ladder are all symmetric Dicke states [3l| . The 
number of motional n p and ionic Hi excitations for each level 
can be inferred from mj, Eq. @. 



To determine the MS states first we rewrite the Hamil- 
tonian (jnj) in terms of the total ionic pscudospin: 



Hj(t) = hg {t) (aJ+ + a f J- ) + HSJ 



(8) 



where J± - EkU and J * = \ E*=i a k ■ Thc MS 
basis consists of the set of the eigenvectors of the two 
commuting operators J 2 and J z , where J 2 = ^(J+J_ + 
J-J+) + J 2 - Therefore each MS state is assigned two 
quantum numbers, respectively j and rrij. Since J 2 com- 
mutes also with the Hamiltonian ([8j), the Hilbcrt space 
factorizes into a set of decoupled chains with different 
values of j. 

The MS states are sketched in Fig. Q]for 6 ions with 
3 excitations. Because the laser pulse couples equally to 
every ion in the trap, the longest chain is comprised of the 
symmetric Dicke states IW 7 ^), wherein a given number 
of ionic excitations n, is shared evenly amongst all thc 
ions in the trap 31(. Each chain is assigned different j 



and is comprised of states with different mj , which varies 
from to j, corresponding to the number of excited ions. 
If rii °f t ne i° ns are excited, then 
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(9) 



Since the total number of excitations is half the number 
of ions, then from Eq. © it follows that for each chain 
mj has a maximum value of 0, i.e. states with nij > are 
inaccessible. By analogy with the traditional angular mo- 
mentum operators, it follows that the number of states in 
a chain is equal to j+1 (j of the states are not accessible), 
from which j can be inferred. The longest chain corre- 
sponds to j = N/2, the next longest to j = N/2 — 1 and 
so on, and overall, there are [JV/2 + 1] chains of different 
length in the factorized coupling scheme. 
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For the following analysis, it will be necessary to go 
further and calculate the couplings in the MS basis. The 
coupling between the neighbors \j, rrij) and \j, rrij — 1) 
follows immediately from the matrix elements of the op- 
erators a* and J_, i.e. 

A™ 3 (t) = g(t)^/n p (j + mj )(j-mj + l) 

= 9(t)y/n p (j + m- N/2)(j + N/2+1- m) . 

(10) 

As illustrated in Fig. [TJ there can be many different 
degenerate MS states with the same values of j and vtij . 
We label the states with a given j ^ N/2 and rrij = 
with \Xj(k)}, where k = 1, . . . , Nj, with Nj = C^, 2 _- — 

^N/2-j-V 



IV. IMPLEMENTATION OF GROVER'S 
ALGORITHM 

The manifold T) 1 which encodes the database elements, 
consists of the states, for which the pseudo-angular mo- 
mentum projection is rrij = 0, Eq. ([5]). This requires that 
the total number of ions TV is even. Half of these are in 
state |1), while the other half are in state |0) (rii = N/2). 
The number of elements TV in the database therefore 
scales exponentially with the number of ions N, Eq. (J7J). 
The elements \x) in Eq. ([2]), which belong to T>, can now 
be written as 

\x)=P x \l 1 ...l N/2 N/2+1 ...0 N ), (11) 

where the subscript x runs over all distinct permutations 
P x of the ions' internal states. The initial state \W) is 
thus a symmetric Dicke state \W^ 2 ) of N ions sharing 
N/2 excitations, 

\W% /2 ) = -7= Yl P * I 11 • • • 1 JV/20jv/2+i ..-On)- 

(12) 

Our proposed experimental procedure consists of four 
simple operations, (i) The ions are firstly initialized in 
the entangled Dicke state (fT2"j). This may be achieved 
using very simple adiabatic passage techniques, involv- 
ing either a pair of chirped laser pulses [3l[ |36| or two 
pairs of delayed but overlapping laser pulses [32J, and 
using global addressing, (ii) Synthesis of the inversion- 
about-average operator is appealingly simple: Myy is a 
single red-sideband off-resonant laser pulse addressing all 
ions in the linear chain, (iii) The oracle query M s is also 
a single red-sideband laser pulse, applied on half of the 
ions. After an appropriate number of iterations (0, the 
system evolves into the marked state |s), which can be 
identified by performing (iv) a fluorescence measurement 
on the entire chain. 



A. Synthesis of the inversion-about-average 
operator 

In most existing proposals for implementing Grover's 
search algorithm using trapped ions the generation of 
Mw requires a large number of concatenated physical in- 
teraction steps, even for moderate register size TV. How- 
ever, by restricting the dynamics to the subspace of the 
overall Hilbert space in which only half of the ions are ex- 
cited, this operator becomes possible to synthesize in only 
a single interaction step. This simplification is achieved 
by taking advantage of the fact that both the Hamil- 
tonian © and the state \W$/ 2 ) are symmetric under 
exchange of any two ions. 

The energies of the MS states do not cross in time so 
that in the limit iS > 1/T the transitions between the 
MS states vanish due to the effect of adiabatic complete 
population return [37} . Each of the MS states acquires 
a phase shift ipj (the index j corresponds to the eigen- 
value of J 2 ) and the unitary propagator within the Dicke 
manifold I? is a product of C^/ 2 -i coupled reflections 

M 

N/2-1 Nj 

j=i k=i 

with Nj = C^ /2 _ 3 - C^ /2 _ J _ l and tp w = ip N / 2 - For a 
given value of the coupling, the detuning may be adjusted 
in order to control the phases ipj. Ideally, these interac- 
tion parameters would be chosen such that ipw ^ (e.g., 

(pw = 7r), while (fj=iN/2 — as this would result in Uw 
being identical to the inversion-about-avcrage operator 
M w (ip w ). 

B. Synthesis of the oracle operator 

The effect of each oracle query M s (cj) s ) is to imprint 
a phase shift <fi s on the marked state \s), whilst leaving 
all other computational basis states unchanged. In gen- 
eral this can be achieved by a multiply-conditional phase 
gate upon the internal state of all N ions in the trap. 
When more than a few ions are involved, this becomes a 
prohibitively complicated operation, which generally re- 
quires an immense number of one- and two-qubit gates 
[13I ] . However, since we do not work in the whole Hilbert 
space, but rather in the manifold V, the oracle operator 
can be implemented in a much simpler fashion - by a sin- 
gle red-sideband laser pulse, addressing uniformly those 
N/2 ions in the trap, which share the excitation of the 
marked state. 

Let's consider an example, when the marked state 
is \s) = 1 111000). Then we address the first three 
ions. Since the initial state \W^ 2 ) and the interactions 
with the laser are symmetric under exchange of the first 
three ions, and the last three, as well, during the exe- 
cution of the algorithm the state of the system is a lin- 
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car combination of the states = \W^j^_ k )\W^ 2 ) , 

k = 0, . . . ,N/2; the marked state is \s) = |$o)- Because 
the Hamiltonian, which describes the oracle call, does not 
change the total number of excitations of the first (and 
the last) three ions, i.e. it does not drive the system be- 
tween l^/jj) and |$fc 2 ) for fci 7^ fc 2 , in the adiabatic limit 
the states |<£>fc) acquire only phase shifts (j>k- Therefore 
the propagator in the manifold T> is formally described 
by the action of N/2 coupled reflections, 



U s = M s 



JV/2-: 

n 

k=l 



(14) 



Hence, as in the case of the inversion-about-average op- 
erator, we need to control N/2 phases. The interaction 
parameters should be chosen such that 4> s = (fiw, while 
4>k=£N/2 = O7 as this would result in U s being identical to 
the oracle operator M s (4> s )- 



C. Phase conditions 

The phase conditions for the propagators Uw and 
U s , derived above, cannot be satisfied exactly, since the 
phases tp^ and 0^ are overly commensurate. However 
one can still perform the algorithm with sufficiently high 
fidelity 



P = l-2|Re(/|A/)|. 



(15) 



Here | /) is the state of the system after the completion 
of the algorithm and |A/) is its deviation due to the 
phases deviation. Fortunately, a Fourier expansion re- 
veals that phase deviation of all odd orders leads to the 
purely imaginary bracket (/|A/) and hence do not con- 
tribute to Eq. (fTS"]) . As a result, the algorithm is less 
sensitive to phase deviations since only those of second 
order have a leading affect on the final state populations. 



D. Numerical Demonstration 
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FIG. 2: (Color online) Simulation of the Grover search algo- 
rithm with N ions and N/2 excitations (register dimension 
C^/ 2 )> f° r N = 6 (upper frame), N = 8 (middle frame) and 
N = 10 (lower frame). The system of ions is assumed to 
be initialized in the symmetric Dicke state \ W^/ 2 )- The laser 

pulses have a Gaussian shape, g(t) = <?oe _t f T . The thin ver- 
tical lines display the timing of the oracle call. The computed 
numerically scaled detunings and peak Rabi frequencies for 
the oracle and the inversion-about-average operator, respec- 
tively, are given in the table below. The marked-state popu- 
lation (the fidelity) of around 99% is obtained after n g = 3, 6 
and 12 steps, respectively, in exact agreement with Grover's 
value (J5]). The numerical simulation includes all off-resonant 
transitions to states with rrij 7^ 0. 
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As a demonstration, we have solved the Schrodinger 
equation numerically for a Gaussian pulse shape and a 
constant detuning 6. Sample results are shown in Fig. 
[2] for N = 6, 8 and 10 ions with 3, 4 and 5 excitations, 
respectively, which imply databases of A/" = 20, 70 and 
252 elements. The fidelity plotted on the vertical axis is 
the time-dependent population of the marked state. The 
system of ions is assumed to be prepared initially in the 
Dicke superposition \ W^ 2 ) of ionic collective states, each 
of which contains exactly N/2 ion qubits in state |1) and 
N/2 in state |0). Each Grover iteration consists of a phase 
shift for the marked state (oracle call), which amounts 
a single red-sideband laser pulse, addressing uniformly 
those N/2 ions in the chain, which share the excitation 
of the marked state, followed by a single red-sideband 



pulse, which addresses uniformly the entire ion chain (the 
inversion-about-average operator). The number of steps, 
for which the algorithm singles out the marked item with 
a probability of around 99%, is respectively n g = 3, 6 and 
12 as predicted by Eq. (0. 



V. CONCLUSIONS 

Despite the intense flurry of theoretical and experimen- 
tal activity in the decade following Grover's original pro- 
posal for a quantum speed-up of unstructured search, a 
large discrepancy still remains between the current exper- 
imental state-of-the-art and what is required for a practi- 
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cally useful quantum search. To highlight this incongru- 
ence, we note that the only physical system in which a 
scalable quantum search has been performed is in a chain 
of trapped ions, and in that experiment, the search-space 
totalled just four elements. An important intermediate 
goal on the long road to performing a practically useful 
quantum search is to demonstrate Grover's algorithm in 
a moderately-sized trapped-ion quantum register. The 
standard approach of building up the many-ion Grover 
operator using a network of single- and two-qubit gates 
is inappropriate for this task, since the required physi- 
cal resources far exceed those available in today's exper- 
iments. By contrast, in this article, we have proposed to 
construct the many-ion Grover operator using only two 
single off-resonant laser pulses, with a suitably chosen 
peak Rabi frequencies and dctunings, which synthesize 
the inversion-about-average and oracle operators, each 
in a single shot. 



A simple recipe for synthesizing the inversion-about- 
average operator was derived by factorizing the over- 
all Hilbcrt space into a series of independent ladders. 
The coupling strengths between the MS states were de- 
termined solely through a consideration of the angular- 
momentum structure of the combined ionic pseudospin. 
The technique proposed in this paper raises the prospect 
of demonstrating Grover's algorithm in a moderately 
sized trapped-ion database comprising up to several hun- 
dred elements, and which scales exponentially with the 
number of ions; this is a necessary step on the path 
to demonstrating a practically useful quantum search, 
which remains a long-term goal. 
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